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Abstract 

Water  management  is  critical  for  Proton  Exchange  Membrane  Fuel  Cells  (PEMFC).  An  appropriate  humidity  condition  not  only  can  improve 
the  performances  and  efficiency  of  the  fuel  cell,  but  can  also  prevent  irreversible  degradation  of  internal  composition  such  as  the  catalyst  or  the 
membrane.  In  this  paper  we  built  the  model  of  water  management  systems  which  consist  of  stack  voltage  model,  water  balance  equation  in  anode  and 
cathode,  and  water  transport  process  in  membrane.  Based  on  this  model,  model  predictive  control  mechanism  was  proposed  by  utilizing  Recurrent 
Neural  Network  (RNN)  optimization.  The  models  and  model  predictive  controller  have  been  implemented  in  the  MATLAB  and  SIMULINK 
environment.  Simulation  results  showed  that  this  approach  can  avoid  fluctuation  of  water  concentration  in  cathode  and  can  extend  the  lifetime  of 
PEM  fuel  cell  stack. 

©  2008  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Proton  Exchange  Membrane  Fuel  Cells  (PEMFC)  (also 
known  as  Polymer  Electrolyte  Membrane  Fuel  Cells)  is 
currently  in  a  relatively  mature  stage  for  ground  vehicle  applica¬ 
tions.  For  PEMFC  the  membrane  must  be  sufficiently  hydrated 
because  its  conductivity  depends  critically  on  the  humidity  level. 
Too  little  water  causes  membrane  drying,  which  increases  the 
ionic  resistance,  and  exacerbates  the  voltage  drop  due  to  ohmic 
losses.  Too  much  water  causes  “flooding”,  which  blocks  porous 
passages  and  thus  reduces  the  transport  rate  of  reactants  to  the 
catalyst  site  [1].  An  inappropriate  humidity  condition  not  only 
can  lower  the  performances  and  efficiency  of  the  fuel  cell  but 
also  can  even  lead  to  an  irreversible  degradation  of  internal  com¬ 
ponent  such  as  the  catalysts  or  the  membrane.  For  examples, 
these  will  lead  to  catalysts  desquamated  and  membrane  dam¬ 
aged  in  the  alternate  wet  and  dry  condition.  Therefore,  water 
management  is  a  critical  issue  for  PEMFC. 
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There  were  many  papers  in  the  literature  discussing  the  fuel 
cell  water  management  issues.  Many  experimental  studies  were 
conducted  to  understand  the  water  transport  phenomena  and  to 
characterize  the  factors  that  affect  the  membrane  water  content 
[2,3].  Shimpalee  et  al.  [4]  studied  a  time  dependence  of  three- 
dimensional  simulation  including  water  phase  change  and  heat 
transfer  of  a  PEMFC  model.  The  overshoot  behavior  has  been 
observed  during  a  change  in  the  electrical  load  during  operation 
with  fixed  flow  rates  of  hydrogen  and  air.  Some  mathematic  mod¬ 
els  have  been  developed  to  optimize  the  fuel  cell  design  in  order 
to  maintain  appropriate  membrane  humidity  of  PEMFC  stack. 
Chen  et  al.  [5]  built  a  one-dimension  model  to  investigate  the 
transient  behavior  of  the  water  transport  across  the  membrane  of 
the  PEMFC.  Alessandro  et  al.  [6]  developed  a  control-oriented 
model  of  the  fuel  cell  system  based  on  static  maps  obtained 
from  a  higher  order  1&1D  model.  Wolfgang  et  al.  [7]  set  up 
a  test  bench  to  measure  the  ohmic  resistance  of  fuel  cell  stack 
and  showed  that  the  humidity  of  membrane  has  a  very  strong 
influence  on  the  electrical  output  of  a  fuel  cell.  Gorgun  et  al. 
[1,8]  presented  a  scheme  to  estimate  membrane  water  content  in 
PEMFC  from  the  measurements  of  voltage,  current,  temperature 
and  reactant  pressure.  Hung  et  al.  [9]  developed  a  first  principle 
one-dimensional  water  and  thermal  management  model  to  gen¬ 
erate  the  I-V  curve.  The  model  considers  the  effects  of  water 
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Nomenclature 

a  water  activity 

Afc  fuel  cell  active  area  (cm2) 

c  water  concentration  (mol  cm-3) 

Dw  diffusion  coefficient 

E0  the  open  circuit  voltage 

F  Faraday’s  number  (C) 

I  stack  current  (A) 

J  rotational  inertia  (kg  m2) 

m  mass  (kg) 

M  molecular  mass  (kg  mol- 1 ) 

n  number  of  cells 

rid  The  electro-osmotic  drag  coefficient 

P  pressure  (Pa) 

R  gas  constant 

R'  electrical  resistance  (Q) 

T  temperature  (K) 

V  voltage  (V) 

V '  volume  (cm3) 

W  mass  flow  rate  (kg  s-1) 

t  time  (s) 

tm  membrane  thickness  (cm) 

v  volume  content  of  dry  air 

y  mole  fraction  of  dry  air 

Greek  letters 
Q  humidity  ratio 

0  relative  humidity 

km  membrane  average  water  content 

am  the  membrane  conductivity 

Subscripts 
act  activation 

a,  air  air 

atm  atmosphere 

H2  hydrogen 

an  anode 

ca  cathode 

cone  concentration 

fc  fuel  cell 

gen  generated 

in  inlet 

inj  inject 

m,  mbr  membrane 
N2  nitrogen 

O2  oxygen 

ohm  ohmic 

out  outlet 

ret  reaction 

rm  return  manifold 

sat  saturation 

sm  supply  manifold 

st  stack 

v  vapor 

w  water 


transport  across  the  membrane,  activation  overpotential,  ohmic 
overpotential,  concentration  overpotential,  pressure  drops,  and 
current  density  distribution  along  the  channel  of  a  PEM  fuel  cell. 

Although  many  researchers  studied  the  mechanism  of  water 
management  and  built  3D,  2D,  ID  or  OD  model,  few  researchers 
explored  how  to  control  the  humidity  of  PEMFC  stack.  Trung 
and  Knobbe  [10]  presented  an  excess  water  removal  method  by 
controlling  each  cell’s  exhaust  so  that  only  one  cell  at  any  given 
time  has  an  open  exhaust  port,  thereby,  ensuring  gas  will  flow 
through  that  cell.  Chen  and  Peng  [11]  built  model  of  humidifier 
system  and  developed  proportional  feedback  control  algorithm 
to  regulate  the  inlet  air  relative  humidity.  A  fuzzy  controller  is 
proposed  for  water  management  of  PEM  based  on  Mamadani 
inference  systems  [12].  In  this  paper,  we  built  model  of  water 
management  system  and  proposed  model  predictive  control  to 
maintain  appropriate  humidity  of  membrane  in  PEMFC.  Model 
Predictive  Control  (MPC)  is  widely  adopted  in  industry  as  an 
effective  means  to  deal  with  large  multivariable  constrained  con¬ 
trol  problems.  In  general,  the  MPC  problem  is  formulated  as 
solving  on-line  a  finite  horizon  open-loop  optimal  control  prob¬ 
lem  subject  to  system  dynamics  and  constraints  involving  state 
and  control  variables.  There  are  some  literatures  discussing  MPC 
methods  of  PEMFC.  Golbert  [13]  developed  MPC  controller 
based  on  the  simplified  nonlinear  model  which  satisfies  load- 
change  demands  and  improves  fuel  efficiency.  Vahidi  et  al.  [14] 
proposed  constraint  MPC  to  prevent  oxygen  starvation  in  a  fuel 
cell  during  rapid  current  transients. 

2.  Model  of  water  management  system 

In  PEMFC  stack  oxygen  and  hydrogen  are  consumed,  and 
water  and  heat  are  generated.  Generally,  a  PEMFC  stack  is  inte¬ 
grated  with  four  auxiliary  systems:  hydrogen  supply  system, 
air  supply  system,  cooling  system  and  humidification  system. 
Water  management  system,  which  is  humidification  system,  is 
used  to  maintain  hydration  of  the  polymer  membrane  and  to 
balance  water  consumption  in  the  system.  The  amount  of  reac¬ 
tant  flow  and  the  water  injected  into  the  anode  and  cathode  flow 
streams  affect  the  humidity  of  the  membrane.  Dry  membranes 
and  flooded  fuel  cells  cause  high  polarization  losses.  As  the  cur¬ 
rent  is  drawn  from  the  fuel  cell,  water  molecules  are  produced 
in  the  cathode  and  dragged  from  the  anode  to  the  cathode  by 
the  hydrogen  protons.  As  the  concentration  of  water  in  the  cath¬ 
ode  increases,  the  concentration  gradient  causes  water  to  diffuse 
from  the  cathode  to  the  anode.  Perturbation  in  fuel  cell  humid¬ 
ity  can  be  caused  by  different  mechanisms:  water  generated 
while  load  increases,  changes  in  the  absolute  and  relative  reac¬ 
tant  pressure  across  the  membrane,  changes  in  air  flow  rate  and 
changes  in  stack  temperature,  which  change  the  vapor  saturation 
pressure  [15].  These  mechanisms  indicate  strong  and  nonlinear 
interaction  among  the  humidity  control  tasks.  A  20-40%  drop 
in  voltage  can  occur  if  there  is  no  proper  humidification  control 
[16]. 

We  build  water  balance  equation  by  investigating  stack  volt¬ 
age  model,  water  distribution  of  membrane  and  water  transport 
process  in  cathode  and  anode.  In  this  paper  direct  injection  of 
water  into  cathode  is  used  for  humidification.  The  following 
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equations  are  stack  voltage  model,  water  distribution  of  mem¬ 
brane  and  water  balance  equation  of  cathode  and  anode. 

2.1.  Water  balance  equation  in  anode 


The  governing  equation  for  hydrogen  and  water  in  the  anode 
can  be  written  as 


dmu2/dt  =  Wn2,in  -  Wu2,out  -  Wu2,rct 

0) 

dmWSLn/dt  =  Wv?an,in  —  fkv,an,out  —  fkv,mbr 

(2) 

The  hydrogen  and  water  flowing  to  the  anode  are 

calculated  by 

lkH2,in  =  (1  —  ^2an,in)lkan,in 

(3) 

fly,  an,  in  =  ^2an,infkan,in 

(4) 

where  £?an,in  is  inlet  flow  humidity  ratio  in  anode. 

The  outlet  hydrogen  and  vapor  mass  flow  rate  is  calculated 
by  the  following  equations. 

fkH2,an,out  =  1/(1  T  ^2an,out)!kan,out 

(5) 

fly, an, out  =  0an,out/(l  +  ^2an,  out)  flan,  out 

(6) 

^2an,out  =  3Tv/3Th2  ’  (wH2  ^H2)/(^v,an^H2) 

(7) 

The  rate  of  water  flow  across  the  membrane  WV,mi br  is  deter¬ 
mined  in  the  membrane  hydrogen  model.  The  rate  of  hydrogen 
consumed  in  the  reaction  ITH2,rct  is  a  function  of  the  stack  current 

WH2,rct  =  MK2(nIst)/(2F)  (8) 

2.2.  Water  balance  equation  in  cathode 

The  mass  continuity  is  used  to  balance  the  mass  of  the  three 
elements:  oxygen,  nitrogen  and  water,  inside  the  cathode  volume 
[17]. 

dniQ2/dt  =  Wo2,in  —  fko2,out  —  fko2,rct  (9) 

dm^2/dt  =  WN2,in  -  l^N2,out  (10) 

dntw  ca/dt  =  fkv,ca,in  —  Wv5ca,out  H”  fkv,gen  fky,mbr  H”  f^inj 

(ID 

where  Wmj  is  the  mass  flow  rate  of  injected  water  from  humidi¬ 
fier.  Inlet  mass  of  oxygen,  nitrogen  and  vapor  is  related  to  mass 
flow  rate  and  humidity  of  inlet  air. 

fko2,in  =  };02/(l  +  ^atm)fkca,in  (12) 

fkN2,in  =  yN2/(l  H-  ^2atm)fkca,in  (13) 

fky,ca,in  =  ^2atm/(l  H-  ^2atm)fkca,in  (14) 

yo2  and  y^2  are  oxygen  mole  fraction  and  nitrogen  mole  fraction 
of  dry  air. 

yo2  =  xo2  •  Mo2/Mfm  (15) 

W2  =(l-xo2)-MN2/Maatm  (16) 


where  xq2  is  volume  content  of  oxygen  in  dry  air.  The  air  molar 


mass  is  calculated  by 

Mfm  =  xq2Mq2  +  (1  —  *o2)MN2 

(17) 

The  humidity  ratio  is  then 

i2atm  =  Mv/Ma  •  VmC7(^atm  -  O 

(18) 

The  pressure  of  cathode  is  the  sum  of  the  partial  pressure  of 
oxygen,  nitrogen  and  vapor. 

Pea  —  P02  +  PN2  +  Pv,ca 

09) 

According  to  ideal  gas  law,  the  partial  pressure  of  oxygen, 
nitrogen  and  vapor  is 

Po2  =  (mo2Ro2Tst)/V^ 

(20) 

Pn2  =  (m^2R^2Tst)/V^ 

(21) 

Pv,ca  =  (^v,ca^y7st)/  Vca 

(22) 

Outlet  gas  flow  rate  is  calculated  by  [15] 

flea,  out  =  ^ca,out(7*ca  —  Rvm,ca) 

(23) 

where  the  mass  flow  rate  of  each  species  out  of  the  cathode  is 
calculated  as 

fko2,out  =  fkca  •  ^02/^ca 

(24) 

fkN2,out  =  fkca  ’  mN2/ mca 

(25) 

fkv, ca, out  =  (Pv,caFcaMv)/(PrstmCa)  ■  fkca 

(26) 

where  mca  =  mo2  +  m^2  +  ( Pv,CaVCaMv)/(RTst )  is  the  total 
mass  of  the  cathode  gas. 

The  rate  of  hydrogen  consumed  in  the  reaction  and  water 
generated  in  the  reaction  is  calculated  by 

Wo2,rct  =  M02(nIst)/(4F) 

(27) 

Wv.gen  =  MH2o(n/st)/(2F) 

(28) 

2.3.  Water  transport  process  in  membrane 

The  water  transport  across  the  membrane  is  achieved  through 
two  distinct  phenomena  [18].  First,  water  molecules  are  dragged 
across  the  membrane  from  anode  to  cathode  by  the  hydrogen 
proton.  This  phenomenon  is  called  electro-osmotic.  Secondly, 
the  gradient  of  water  concentration  across  the  membrane  results 
in  “back-diffusion”  of  water,  usually  from  cathode  to  anode. 
The  water  concentration  is  assumed  to  change  linearly  over 
the  membrane  thickness.  Combining  the  two  water  transport 
mechanisms,  the  water  flow  across  the  membrane  from  anode 
to  cathode  [15]  is 

lFv?mem=3Tv^lfc^(^d7st/ F  ~  Dw  '  (cv,an  —  cv,an)/ hn)  (29) 

The  electro-osmotic  drag  coefficient  n d  and  the  diffusion 
coefficient  Dw  vary  with  water  content  in  the  membrane  [19,20] . 

nd  =  0.0029A.2  +  0.05V  -  3.4  x  10"19  (30) 

Dw  =  Dx  exp(24 16(1/303  -  1/V))  (31) 
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where 

'  10“6,  v  <  2 

„  10“6(1  +  2(V  -  2)),  2  <  v  <  3 

D\  —  ,  (3 

10_6(3  —  1.67(A.m  —  3)),  3  <  Xm  <  4.5 

1.25  x  10“6,  km  >  4.5 

Membrane  average  water  content  Xm  is  calculated  by  [19,20] 


— 

where 


0.043  +  17.81<2m—  39.85^  +  36.0 a^,  0  <  am  <  1 
14  +  lA(am  —  1),  1  <  am  <  3 

(33) 


am  —  (fl an  H”  aCSL)/2  (34) 

Water  activities  in  anode  and  cathode  <2an  and  aCSL  is  calculated 
by 


^an  —  ^v,an/ ^sat,an>  ^ca  —  ^v,ca  /  ^sat,c 


(35) 


2.4.  Stack  voltage  model 


Va  =  (-1.618  x  10"57’+  1.618  x  10“2)( - — - 

V0.1173  +  psat 

x(1.8  X  10“4r  -  0.166)  ( - — - ) 

VO-1173  +  psat/ 

+(-5.8  X  10“4r  +  0.5736)  (41) 

Ohmic  voltage  loss  Vohm  is  calculated  as  follows 

Vohm  =  ^'ohm  =  —  (42) 

where  tm  is  the  thickness  of  membrane,  am  is  the  membrane 
conductivity  and  calculated  as  follows 

<%  =  (.,  exp  (43) 

in  which  Z?2  is  constant  [21],  Z?i  is  function  of  water  content  of 
membrane 

bi  =  0.005139V  -  0.00326  (44) 

Concentration  voltage  loss  Vconc  is  calculated  as  follows 


The  stack  voltage  is  calculated  as  a  function  of  stack  current, 
cathode  pressure,  reactant  partial  pressures,  stack  temperature 
and  membrane  humidity.  There  are  a  lot  of  stack  voltage  model 
in  the  literature.  We  adopt  the  semi-empirical  model  in  this 
paper  [21]. 


Tconc  —  *(C2*/imax)C3  (45) 

where  ci,  c  \  and  /max  are  constant  [21]. 

3.  MPC  based  on  RNN 


vfc  =  Eo-  Ka  -  V'hm  ~  V'onc  (36) 

where  E$  is  the  open  circuit  voltage,  Vact  is  the  activation  voltage 
loss  resulted  from  the  need  to  move  electrons  and  to  break  and 
form  chemical  bonds  in  the  anode  and  cathode,  V0hm  is  the  ohmic 
voltage  loss  resulted  from  the  resistance  of  the  polymer  mem¬ 
brane  to  the  protons  transfer  and  the  resistance  of  the  electrode 
and  the  collector  plate  to  the  electrons  transfer,  VCOnc  is  the  con¬ 
centration  voltage  loss  resulted  from  the  drop  in  concentration 
of  the  reactants  as  they  are  consumed  in  the  reaction. 

Eo  is  the  function  of  stack  temperature  T,  pressure  of  hydro¬ 
gen  Ph2  and  pressure  of  oxygen  Po2  [15,21] 

E0  =  1.229  -  0.85  x  10 ~3(T  -  298.15)  +  4.3085 

x  10-5r  [ln(pu2)  +  0.5 \n(po2)\  (37) 


MPC  is  widely  adopted  in  industry  as  an  effective  means 
to  deal  with  large  multivariable  constrained  control  problems. 
At  each  sample  instance  MPC  problem  has  a  quadratic  objec¬ 
tive  function  and  subject  to  system  dynamics  and  constraints 
involving  state  and  control  variables.  So  these  form  the  Quadratic 
programming  (QP)  problem  at  each  sample  instance.  As  the  QP 
optimization  problem  depends  on  the  current  state,  the  imple¬ 
mentation  of  MPC  requires  the  online  solution  of  a  QP  at  each 
sample  instance.  High  computational  burden  in  solving  QP  prob¬ 
lem  is  a  major  obstacle  when  we  apply  MPC  to  industrial 
process.  To  solve  this  problem,  a  novel  MPC  based  on  RNN 
is  proposed. 

3.1.  Formulation  of  MPC 


where 


Ph2  = 


(mn2  Ru2Tst) 
~V*n 


Activation  voltage  loss  Vact  is  calculated  as  follows 
Vact  =  Vo  +  Va(l  -  e-C1') 


(38) 


(39) 


where  ci  is  parameter  [15,21],  Vo  and  Va  can  be  calculated  as 
follows 

V0  =  0.279  -  8.5  x  10“4(7’  -  298.15)  +  4.308  x  10“57’  • 


In  |  Pca  M+lln 
1.01325  J  2 


0.01 173(/?ca  Psat)\ 
1.01325  ) 


(40) 


Consider  the  following  receding  horizon  formulation: 

N 

min  =  VOV  +  k)~  yr)T  Q(y(t  +  k)~  vr) 

U  Z ' 

k=  1 

+{u(t  +  k)  —  uv)T  R(u(t  +  k)  —  Mr)  (46) 

s.t.  x{t  T  k  1)  =  Ax(t  T  k)  -T  Bu(t  T  E) 
y(t  4-  k  -f  1)  =  Cx(t  T  E)  T  Du(t  -f  E) 
umm  —  u(t  +  E)  <  Mmax 

where  x{t  +  E)  e  Rn  is  the  expected  value  of  the  state  vari¬ 
able  at  the  time  t  +  k  as  predicted  at  the  time  t ,  y(t  +  k)eRp 
is  the  corresponding  output  variable,  u(t  +  E)  e  Rm  is  manip¬ 
ulated  variables  considered  in  the  optimization  program  at 
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the  time  t  +  k  as  predicted  at  the  time  t.  N  is  predictive 
horizon,  yr  and  ur  are  reference  trajectory  of  output  vari¬ 
able  and  manipulated  variable,  respectively.  QeRpxp  and 
R  e  Rmxm  are  weight  matrices.  x(t  +  k  +  1)  =Ax(t  +  k)  +  Bu(t  +  k) 
and  y(t  +  k  +  1 )  =  Cx(t  +  k)  +  Du(t  +  k)  are  plant  model. 

It  is  well  known  that  MPC  solves  this  problem  at  each  sample 
instance  t  by  computing  all  next  N  values  of  the  manipulated 
variables  u(t+k).  According  to  system  model,  we  can  get 

k 

x(t  +  k  +  1)  =  Akx(t )  +  ^^Ak~l  Bu(t  +  i)  y(t  +  k  +  1) 
i=  1 

k—  1 

=  CAk~xx(t)  +  Buit  +  i)  +  Du(N) 

i= 1 

(47) 

Substituting  Eq.  (47)  into  optimization  problem  (46),  we  have 
the  form  of  a  QP  problem 

min  | -UT HU  +  ( x{t)T S  +  T)U  1 
t/  \2  J  (48) 

S.t.  t/min  ^  ^  t^max 

where  x(t)  is  value  of  state  variable  at  time  t,  U-  [u(t),u(t  +  1), 
...  ,u(t  +  N)]T  is  manipulated  variable  set  in  next  N  step, 
HeRmNxmN ?  S  e  RnxmN ,  T  e  RmN ,  e  e  RmNxmN  is  fixed  value 
matrix  and  can  be  easily  obtained  from  the  plant  and  Q ,  7?  (for 
details,  see  [22]). 

3.2.  7?A7V  optimization 


It  can  be  seen  that  the  state  equation  described  by  Eq.  (49)  can  be 
implemented  by  an  analog  circuit  with  a  single-layer  structure. 
The  proposed  MPC  based  on  RNN  method  can  be  summarized 
as  in  the  following: 

(1)  Initialization.  Let  k-  1.  Set  terminal  time  rend,  the  value 
of  state  variable  x(k)  in  plant  model,  predictive  horizon 
N,  control  horizon  Nc ,  sample  time  t ,  weight  matrix  Q 
and  P,  computation  time  At  of  RNN  at  each  sample  time. 
According  to  the  plant  model,  calculate  matrix  H ,  S ,  T  and 
W=R  =  H~ 1 . 

(2)  Calculate  RNN  parameters.  Calculate  f=S-x(k)  +  T, 
q  =  a=  —  H~xf 

(3)  RNN  dynamic  optimization.  According  to  W,  R ,  q  and  a , 
RNN  start  to  run,  after  time  At,  we  get  stable  value  U*. 

(4)  Prediction  of  plant  model.  Get  first  m  value  in  if  as  the 
value  of  manipulated  variable  and  calculation  next  step 
value  x(t+  1)  of  state  variable. 

(5)  If  k<T& nd,  set  k=k+  1,  return  step  2,  otherwise,  program 
end. 

To  study  optimization  capability  of  the  proposed  MPC 
based  on  RNN,  convergence  property  of  RNN  at  each  sample 
instance  must  be  investigated.  According  to  theorem  in  [23], 
our  main  result  on  the  proposed  MPC  based  on  RNN  is  the 
following. 

Theorem  1.  The  state  trajectory  of  the  proposed  recurrent  neu¬ 
ral  network  is  exponentially  convergent  to  a  solution  of  Problem 
(48). 


At  each  step  x(t)  is  already  known  in  QP  problem  (48),  so  we 
apply  RNN  proposed  by  Xia  et  al.  [23]  to  solve  this  QP  problem. 
Consider  the  Lagrangian  function  of  problem  (48) 

L(x,  y,  rj)  =  ^xtHx  +  ((x(t)T S  +  T)U  -  kr(EU  -  rj) 

where  X  e  RmN  is  referred  to  as  the  Lagrange  multiplier  and 
t]  g  X=  {X  e  RmN\um[n  <X<  um  in}.  According  to  Saddle  point 
theorem,  we  assume  that  the  value  of  X  and  r\  are  X*  and  rf 
when  v*  is  an  optimal  solution  of  the  problem  (48). 

In  [23],  an  approach  to  transform  QP  problem  into  RNN 
optimization,  whose  dynamic  equation  is  defined  as  follows. 
State  equation 

—  =  X{Px(Wu  +q-u)-Wu-q}  (49) 

d t 

Output  equation 


U(t)  =  Rv(t)  +  a  (50) 

where  W=  R  =  Q~l  is  matrix,  q  =  a  =  Q~l(x(l)T R  —  T)  is  vector, 
v  e  RmN  is  the  state  variable  of  RNN,  X  >  0  is  a  scaling  constant. 

Px(v)  =  [Px(t>i),  Px(v  2),  •  •  • ,  Px(vmN)]T  and  for;  = 

1,  . . . ,  m  *  N 


Px(vt)  = 


wmin  (0  Vi  <  wmin  (0 

Vi  Vf  <  umSLX(i) 

wmax(0  Vi  >  UmSLX(i) 


(51) 


Theorem  2.  Assume  that  W  =  EQ~JET  and  q  =  EQ~\x{t  +  1)T 
R  —  T),  the  output  trajectory  of  the  proposed  RNN  converges 
globally  to  a  unique  optimal  solution  of  Problem  (48)  at  each 
step  within  a  finite  time.  Moreover,  it  has  the  following  converge 
rate 

\\x(t)  -  x*\\2  < -f  —  Vt  >  to 
X(t  -  to) 

where  y  is  a  positive  constant. 

4.  Discussion  and  results 

MPC  based  on  RNN  methods  are  used  to  improve  water  man¬ 
agement  in  PEMFC.  The  control  objective  of  water  management 
is  to  keep  water  concentration  of  cathode  in  a  fixed  value  and 
avoid  fluctuation  of  water  concentration  in  cathode. 

According  to  the  nonlinear  model  described  in  Sec¬ 
tion  2,  the  state  equations  involve  five  state  variables  v  = 
[mu2,  mw?an,  mo2,  m^2,  mw?ca]r.  The  parameters  used  in  this 
model  are  given  in  Table  1.  Most  of  the  parameters  are  based 
on  the  75  kW  stacks  used  in  the  FORD  P2000  fuel  cell  proto¬ 
type  vehicle  [15,24].  This  stack  consists  of  381  single  fuel  cells, 
the  membrane  of  single  fuel  cell  is  Nation  117.  All  nonlinear 
models  are  implemented  in  Matlab/Simulink  environment.  The 
Simulink  block  is  shown  in  Fig.  1 .  This  model  is  capable  of  cap¬ 
turing  the  effects  of  transient  humidity  on  the  fuel  cell  voltage. 
Even  the  model  has  not  been  validated  with  an  actual  experi- 


Table  1 

Parameter  used  in  the  model 


L.  Zhang  et  al.  /  Journal  of  Power  Sources  180  (2008)  322-329 


327 


Fig.  2.  Simulation  result  without  control. 


Symbol  Variable  Value 


Mh2 

Hydrogen  molar  mass 

2.016x10  3  kg  mol' 

MQj 

Oxygen  molar  mass 

32  x  10_3kgmol_1 

Mn2 

Nitrogen  molar  mass 

28  x  10-3  kg  mol-1 

*h2 

Hydrogen  gas  constant 

4124.3  J  (kg  K)”1 

*n2 

Nitrogen  gas  constant 

296.8  J  (kg  K)"1 

Ro2 

Oxygen  gas  constant 

259.8  J  (kg  K)-1 

Ry 

Vapor  gas  constant 

461.5  J  (kg  K)-1 

Ra 

Air  gas  constant 

286.9  J  (kg  K)-1 

P atm 

Atmospheric  pressure 

101.325  kPa 

n 

Number  of  cells  in  fuel  cell  stack 

381 

Ac 

Fuel  cell  active  area 

280  cm2 

Fan 

Anode  volume 

0.005  m3 

Vca 

Cathode  volume 

0.01  m3 

F 

Faraday  number 

96485 

£2  an, in 

Humidity  ratio  of  hydrogen 

0 

0atm 

Relative  humidity  of  air 

0.5 

mental  system,  the  model  predicts  transient  behavior  similar  to 
that  reported  in  the  literature  [15]. 

Firstly,  we  see  the  simulation  result  when  we  do  not  control 
mass  flow  rate  of  injected  water  from  humidifier  W\nj .  During 
0-3  s  time  period  stack  current  is  200  A  and  steps  up  to  300  A 
during  3-6  s  time  period.  Mass  flow  rate  of  injected  water  is 
0.06  kg  s-1  during  0-6  s  time  period.  Simulation  result  is  shown 
in  Fig.  2.  From  Fig.  2  we  see  that  water  mass  in  cathode  is 


about  0.0021  kg  during  0-3  s  time  period  and  drop  to  0.001 86  kg 
during  3-6  s  time  period,  stack  voltage  is  about  264  V  during 
0-3  s  time  period  and  about  208  V  during  3-6  s  time  period.  We 
assumed  that  the  inlet  airflow  rate  is  well  controlled  in  all  cases. 
The  inlet  airflow  rate  is  0.4  kg  s_1  during  0-3  s  time  period  and 
0.6  kg  s-1  during  3-6  s  time  period.  If  we  do  not  control  mass 
flow  rate  of  injected  water  when  stack  current  step  up  from  200 
to  300  A,  water  humidity  in  cathode  cannot  maintain  a  stable 
value  and  stack  voltage  decreased  with  a  large  amount. 

According  to  the  nonlinear  model  shown  in  Fig.  1,  we  adopt 
PID  controller  to  control  water  mass  in  cathode.  Control  vari- 


Fig.  1.  Simulink  block  of  water  management  model. 
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Fig.  3.  Simulation  result  of  PID  controller. 


Fig.  5.  Simulation  result  of  MPC. 


able  is  mass  flow  rate  of  injected  water  from  humidifier  W-my 
Because  it  is  hard  to  measure  water  mass  in  cathode,  we  designed 
an  observer  of  water  mass  in  cathode  according  to  water  man¬ 
agement  model  similar  as  [25].  During  0-3  s  time  period  stack 
current  is  200  A  and  steps  up  to  300  A  during  3-6  s  time  period. 
Simulation  result  is  shown  in  Fig.  3 .  We  notice  that  water  mass  in 
cathode  can  be  kept  at  a  stable  level  about  0.0021  kg,  stack  volt¬ 
age  is  about  264  V  during  0-3  s  time  period  and  about  240.5  V 
during  3-6  s  time  period,  and  the  response  time  is  about  0.5  s 
when  stack  current  steps  up. 

When  we  linearize  this  system  at  the  nominal  operating 
point  with  stack  current  as  200  A  and  stack  temperature  as 

x  =  Ax  +  Bu 

338  K,  we  can  get  state  space  equation  ,  where 

y  =  Cx  +  Du 

x  ~  [mH2,  mw,an,  WIO2,  mN2,  Wiw>ca]r*  y  is  the  water  mass  in 
cathode  and  u  is  the  flow  rate  of  injected  water  from  humidifier. 
The  reference  output  is  Yv  =  0.0021  kg,  P  =  1,  Q  =  10,  predictive 
horizon  N=  3,  sample  time  t  =  0.1  s.  We  used  model  predictive 
control  based  on  RNN.  In  this  controller  control  variable  is  the 


time  (s) 


flow  rate  of  injected  water  in  humidifier  and  output  variable  is 
water  mass  in  cathode.  The  RNN  has  3  state  variables.  Tran¬ 
sient  behaviors  of  RNN  state  at  each  sample  instance  is  shown 
in  Fig.  4.  From  Fig.  4  we  see  that  x  is  the  state  variable  of  RNN  at 
each  sample  instance  and  convergence  time  is  about  0.005  s  and 
20  times  shorter  than  sample  time.  So  the  computation  time  of 
RNN  optimization  is  enough  to  solve  the  MPC  problem  at  each 
sample  time.  Simulation  result  of  MPC  controller  is  shown  in 
Fig.  5.  From  Fig.  5  we  see  that  water  mass  in  cathode  can  keep 
the  fixed  value  about  0.0021  kg,  stack  voltage  is  about  264  V 
during  0-3  s  time  period  and  about  240.5  V  during  3-6  s  time 
period,  and  the  response  time  is  about  0.15  s  when  stack  current 
steps  up.  When  the  water  mass  in  cathode  is  0.0021  kg,  humidity 
in  cathode  is  about  90-100%.  Moreover,  we  simulate  contin¬ 
uous  change  of  the  multiple  stack  currents.  During  0-3  s  time 
period  stack  current  is  100  A,  steps  up  to  150  A  during  3-6  s  time 
period,  steps  up  to  200  A  during  6-9  s  time  period,  steps  down 
to  150  A  during  9-12  s  time  period,  steps  up  to  250  A  during 
12-15  s  time  period  and  steps  up  to  300  A  during  15-18  s  time 
period.  Simulation  result  is  shown  in  Fig.  6.  So  MPC  controller 
can  keep  humidity  of  cathode  in  fixed  value.  Moreover,  MPC 
controller  has  shorter  response  time  than  PID  controller.  From 
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Fig.  4.  Transient  behaviors  of  RNN  state  at  each  sample  instance. 


Fig.  6.  Simulation  result  of  MPC  in  different  stack  current. 
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the  comparisons  of  three  cases:  without  control,  PID  control  and 
MPC  we  can  conclude  that  MPC  controller  can  avoid  fluctuation 
of  humidity  in  cathode  and  prevent  irreversible  degradation  of 
internal  component  such  as  the  catalyst  or  the  membrane. 

5.  Conclusion 

In  this  paper  a  detailed  model  of  water  management  in 
PEM  fuel  cell  is  built  and  its  simulation  is  implemented  in 
the  Matlab/Simulink  environment.  Based  on  these  model  and 
simulation,  we  proposed  MPC  approach  to  control  water  con¬ 
centration  in  cathode  and  reduce  fluctuation  of  humidity  of 
cathode.  MPC  optimization  at  each  sample  instance  is  achieved 
using  RNN  optimization.  Simulation  results  showed  that  MPC 
controller  based  on  RNN  has  shorter  response  time  than  PID 
controller  and  can  avoid  fluctuation  of  water  concentration  in 
cathode.  So  this  method  can  be  applied  to  real-time  control  in 
water  management  of  PEMFC. 
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